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Abstract 



We calculate the D-wave heavy quarkonium production at fixed target ex- 
periments under the NRQCD factorization formalism. We find that the color 
octet contributions are two orders of magnitude larger than color-singlet con- 
tributions if color-octet matrix elements are taken according to the NRQCD 
velocity scaling rules. Within the theoretical uncertainties, the prediction for 
the production rate of 2 D-wave charmonium state agrees with the prelim- 
inary result of E705 and other experiments. Searching for the 1 D-wave 
state ^p{3770) is further suggested. 
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Typeset using REVT[TiX 



Studies of heavy quarkonium production in high energy colhsions provide important in- 
formation on both perturbative and nonperturbative QCD. Recent progress in the this area 
was stimulated by the experiment resuhs of CDF at the Fermilab Tevatron. In the 1992-1993 
run, the CDF data |]1| for the prompt production rates of ip and ip' at large transverse mo- 
mentum were observed to be orders of magnitude larger than the lowest order perturbative 
calculation based on the color-singlet model. At the same time, a new framework in treating 
quarkonium production and decays has been advocated by Bodwin, Braaten and Lepage in 
the context of nonrelativistic quantum chromodynamics (NRQCD) In this approach, the 
production process is factorized into short and long distance parts, while the latter is associ- 
ated with the nonperturbative matrix elements of four-fermion operators. This factorization 
formalism gives rise to a new production mechanism called the color-octet mechanism, in 
which the heavy-quark and antiquark pair is produced at short distance in a color-octet 
configuration and subsequently evolves nonperturbatively into physical quarkonium state. 
The color-octet terms in the gluon fragmentation to J/tplip') have been considered to ex- 
plain the J/ip{ip') surplus problems discovered by CDF [^Q. Taking the nonperturbative 
(O'l^'^ i^Si)^ and (^5*1)^ as input parameters, the CDF surplus problem for J and ip' 
can be explained as the contributions of color-octet terms due to gluon fragmentation. In the 
past few years, applications of the NRQCD factorization formalism to J/ip{ip') production 
at various experimental facilities have been studied 0. 

Even though the color-octet mechanism has gained some successes in describing the 
production and decays of heavy quark bound systems, it still needs more effort to go be- 
fore finally setting its position and role in heavy quarkonium physics. (For instance, the 
photoproduction data from HERA |^ puts a question about the universality of the color- 
octet matrix elements 0, in which the fitted values of the matrix elements {O'l^^ {^Sq)) and 
{pl^^ i^Pj)) are one order of magnitude smaller than those determined from the Tevatron 
data So, we must find other processes to test the color-octet mechanism in the heavy 
quarkonium production. 

In our previous studies [p|-[T2|l, we propose the D-wave heavy quarkonia production as 



a crucial test for the color-octet mechanism. We have calculated the D-wave quarkonium 
production via gluon fragmentation at the hadron colliders ||10[ , in the decays [O , and the 



Z^-wave charmonia production in B decays [0. All these results show that the color-octet 
mechanism is crucially important to D-wave charmonium production, and the color-octet 



^Possible solutions for this problem have been suggested recently in 
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contributions are found to be over two orders of magnitude larger than the color-singlet 
contributions. 

In this paper, we analyze the D-wave heavy quarkonium production in fixed target exper- 
iments, at energies in which the color-octet gluon fragmentation is not dominant. However, 
the analysis of S-wave charmonium and bottomonium production in fixed target experiments 
show that the color-octet contribution is important in removing large discrepancies between 
experiments and predictions from the color-singlet model |T^. We can further expect that 



the color-octet contribution may be dominant in the production of D-wave quarkonium pro- 
duction in fixed target, because within the color-singlet model the D-wave heavy quarkonia 
production is always suppressed by the second derivative of the wave function at the origin. 
But this suppression can be avoided in the color-octet model. Therefore, the D-wave heavy 
quarkonium production may provide another test for the color-octet mechanism. In addi- 
tion, on the experimental side, there are some clues for the D-wave 2 charmonium state 
in E705 300 GeV vr^- and proton-Li interaction experiment In this experiment there 
is an abnormal phenomenon that in the J/ipTi^TT' mass spectrum, two peaks at 'i/'(3686) 
and at 3.836 GeV (given to be the 2 state) are observed and they have almost the same 
height. Obviously, this situation is difficult to explain based upon the color-singlet model. 
However, it might be explained within the scope of NRQCD. 
In NRQCD the Fock state expansion for ^Dj states is 

fDj) = Oil)\QQi'Dj, D) + Oiv)\QQCP,r, 8)g) + Oiv')\QQi'Su 8 or l)gg) + ■■■. (1) 

In the above expansion, the contributions to the NRQCD matrix elements for the production 
of the D-wave charmonium from the three terms are the same order of v according to the 
NRQCD velocity scaling rules. So, in the factorization formula of D-wave quarkonium 
production, all of these three terms should be considered. Accordingly, the production cross 
section for the physical D-wave quarkonium state 6j in hadron process 

A + B^6j + X (2) 

can be written as 

1 

a{A + B ^ 5j + X) = J2 J dxidx2fi/A{xi)fj/B{x2)cr{ij 5j), (3) 

i,j 



^Here, the symbol 5j denotes the physical spin-triplet D-wave heavy quarkonium state. The 
notation "^^^^Dj represents the cc or bb pair configurations with angular momentum L = 2. 
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cT{tj^5j) = Y.F{zj^n){0'^'). (4) 

n 

Here, n denotes cc pair configuration in the three Fock states of Eq. (P (including angular 
momentum "^^^^Lj and color index 1 or 8). Fiij — > n) is the short distance coefficient for 
the subprocess ij n. (O^^^ is the long distance non-perturbative matrix element which 
represents the property of the cc pair in n configuration evolving into the physical state 5j. 
The short distance coefficient F can be calculated by using perturbative QCD in expansion 
of powers of a^. The long distance matrix elements are still not available from the first 
principle at present, which can are obtained by fitting the theoretical prediction with the 
experimental result in literatures. By the naive NRQCD velocity scaling rules, these three 
matrix elements are of the same order in powers of v. Their scaling properties are, 

0{^{^Dj))r^Mlv\ (0t-'i'S,))^My, M'{^P'j))^My. (5) 



That is to say that the color-octet and color-singlet processes are of the same order in 
powers of v. However, the color-octet processes are enhanced by an order of over the 
color-singlet processes. So, in the hadropro duct ion of D-wave heavy quarkonium at fixed 
target experiments the color-octet contributions should be included. 

First, we will calculate the color-singlet contribution to ^Dj quarkonium hadroproduc- 
tion. The leading order color-singlet contribution comes from the gluon-gluon fusion process 
gg — > ^Djg. The production rate of this process can be calculated by making use of the 
covariant formalism. The general form of the wave function of the spin-triplet heavy quarko- 
nium state (with angular momentum ^Lj) can be written as 

where e'-''^ is the polarization vector associated with the spin-triplet states. A^(pi) and 
A^(p^) are positive energy projection operators of quark and antiquark . 

A+(pi) = , A_(p2) = . (7) 

For the D-wave function, $(P, g) must be expanded to the second order in the relative 
momentum q. The first and second derivatives of the wave function are 

^a{q) = ^^J2iJM\lsLm)[M,^a4^^\M+ /') + M,/W(M+ /')7J^L„^(g); (8) 
^apiq) = ^^T.iJM\lsLm)^a /'\M+ r)lp^Lm{q). (9) 
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After integrating over the relative momentum g*, the orbit angular momentum part of the 
wave function will depend on the radial wave function or its derivatives at the origin Rs{0), 
R'p{0), and -R'^(O), respectively, for S-wave (L = 0), P-wave (L = 1), and D-wave (L = 2) 
states, 

'^'^-^oo(q) = ^i?5(0), (10) 



'^'^g"V^l™(q) = ^e^&^(0), (11) 



'^'^g"g^^2™(q) = el\^R'l,{0), (12) 



(27r)3 

where the polarization tensor's label m is magnetic quantum number. For the spin-triplet 
case where J = 1, 2, 3, using explicit Clebsch-Gordan coefficients, there are following 



relations for the three cases [|T5|. 



-^9.. - ' (13) 

g{2J.|ls2m)6!;'''6('> = -^(4-J)e,^,,,p-s-' + 4i->£™,,,p-9"'), (14) 

Y.{3J,\ls2m}it';'i';>^Af^- (15) 

Here, e^, ea/3, ^app are the spin-one, spin-two and spin-three polarization tensors which obey 
the projection relations [|T^ 

E4'"H"^^ = (-^7.^ + |g)-^./3, (16) 



m 



1 1 

^«™/3' = l^^Paa'Vpp' + Va(3'Vpa'] " -Va/sVa'P', (17) 



m 



1 



m 



+ 'Pal3''Pl3p''Ppa' + 'Pap''Pf3l3''Ppa' + 'Pap''Pl3a''Ppl3'] 
— i'Pal3'Ppa''P/3'p' + 'Pal3'Ppl3''Pa' p' + a({P pp'^P a' 

+ VapVpa'ViS'pi + VapVppiVa'p' + VapV ppiVa' fS' 
+ 'Pf3p'Paa''Pf3'p' + 'P^p'Pal3''Pa' p' + 'Pi3p'Pap''Pa' 13'] 
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So, the differential cross section for the color-singlet process gg — > 5jX can be calculated 
as an expansion of the form in Eq.(^, 

da{gg-^6jg), _ p. _ 3 n /.o^/^s 



Figg^-'Djg)x{Ol'i-'Dj)), (19) 

where the color-singlet matrix element (Oi'i^Dj)^ can be related to the second derivative 
of the nonrelativistic radial wave function at the origin |i?^(0)p for D-wave by 

,5,,,3n A\ _ 15(2J + l)iV, 



OVi'Dj)) = -'^'^ ' ^ |i?^(0)|^. (20) 

The calculation of the short distance coefficient Fi^gg ^Dj)g is straightforward. Because 
the results are lengthy, we give the expressions in the Appendix. As in the case of gluon frag- 
mentation to color-singlet ^Dj [|T^, gg — ^ ^Djg processes also have the infrared divergences 
involved, which are associated with the soft gluon in the final state. In our calculations of 
the cross section for these processes, we follow the way in Ref. |T0| by imposing a lower cutoff 
A on the energy of the outgoing gluon in the quarkonium rest frame. The cutoff A can be 
set to be mq to avoid large logarithms in the divergence terms. 

For the color-octet contributions in the gg, gq, qq subprocesses, we readily have @,|13 



aigg ^ - t|3^^(-.- - (oi-^CP.)) + 4 {o^-^eP.))l (21) 

^igq ^ Sj)s = 0, (22) 
- Sj)s = ^^(^^^(^1^2 - 4mys) (o's^CS,)) . (23) 

Here ^/s is the center-of-mass energy, and is normalized at the scale 2mQ. 

Our numerical results are displayed in Fig.l to Fig. 4. We use the Gliick-Reya-Vogt (GRV) 



LO parton distributions for the proton and the pion [^,|T^. We set the renormalization scale 
to be 2mQ. 

In Fig.l and Fig.2, we plot the total cross section of D-wave heavy quarkonium for 
proton-nucleon collisions for xp > 0, where we choose the color-octet matrix elements to be 
related to the color-singlet matrix elements according to the naive NRQCD velocity scaling 
rules, 



olh'Si)) (ofeP,)) (ofi^D, 



M| M5 MJ 

0lh'S^)) (ofi'P,)) lo'}{^D,] 



Ml Ml Ml 



(24) 
(25) 
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And the heavy quark spin symmetry relations for the color-octet matrix elements have been 
used, 



O's^CS,)) (Orr^i)), (26) 



Ol'CPi)) ^^^{O^'CPi)). (27) 

The values of color-singlet matrix elements are obtained from the potential model calculation, 
\R'[j{0)\l = OMbGeV^, |i?'/,(0)|f = OmTGeV^ P]. Fig.l is the cross section for the D-wave 



charmonium and Fig. 2 for the D-wave bottomonium respectively. From these figures, we 
can see that the color-octet contributions are two orders of magnitude larger than the color- 
singlet contributions. 

It must be noted that all of the cross sections are expressed in terms of heavy quark mass 
as rric and rrif,, for which we take values as 

m, = l.5GeV, = A.^GeV. 

If we choose the charm quark mass as one half of the charmonium mass (e.g., 5^ ~ 3.8GeV^), 
the theoretical prediction of the cross section will be reduced by an order of magnitude for 
the color-singlet contributions and by a factor of five for the color-octet contributions at a 
typical energy scale ^/s = 24GeV. Compared with this large uncertainty, the variation due 
to the choice of parton distribution functions and as{fi) is neghgible. 

In Fig. 3 and Fig. 4, we plot the cross sections of D-wave charmonium and bottomonium 
states in pion-nucleon collisions for xp > 0. Similar to the case of proton-nucleon collisions, 
in pion-nucleon collisions we can see from these figures that the color-octet contributions are 
two order of magnitude larger than the color-singlet contributions. 

Among the three triplet states of D-wave charmonium, 62 is the most promising candidate 
to discover firstly. Its mass falls in the range of 3.810 ~ 3.840 GeV in the potential model 
calculation, that is below the DD* threshold but above the DD threshold. However, the 
parity conservation forbids it decaying into DD. It, therefore, is a narrow resonance. Its 
main decay modes are expected to be, 

^c^j/^TTTT, --'PMJ = 1,2), 5l^?>g. (28) 



The hadronic transition rate of 82 — > J/^'t^^t^ is estimated to be |jTT| 

T{51 J/^/'TT+TT") = V{^Di J/^/'TT+TT") ^ 46 kcV. (29) 

For the El transition 5l —>■ '^Pj'~f{J = 1, 2), using the potential model with relativistic effects 
being considered IT^, we find 



7 



r(52^ Xcil) = 250 keV, V{51 ^ Xc2l) = 60 keV, (30) 
where the mass of 82 is set to be 3.84GeV. As for the 82 — ^ 3(7 annihilation decay, an estimate 



gives [jT5| 

V{51 ^ 3^) = 12 fceK 
From (|2|), dg), and (|1]), we find 

rto,(52^) ^ r(52^ J/^vrvr) + V{51 ^ Xci7) + r(52' ^ Xc27) + r(52^ ^ 3<7) 



(31) 



390 /ceK 



and 



B{6'2 J/^TT + TT") 0.12. 



(32) 



(33) 



Comparing ([33| ) with -B('?/'' ^ J/ipTr^Ti^) = 0.324 ± 0.026, the branching ratio of 62 
J/ip7T^7T~ is only smaller by a factor of 3, and therefore the decay mode of 62 — > J/iIjtt'^tt~ 
could be observable in various experiments, such as at hadron colliders and at fixed target 
experiments. 

Multiplied by the assumed branching ratio of S2 — > J/ipn^Ti^ above, we can estimate 
the production rate of 62 in pion-nucleon collisions. From Fig. 3, at ^ 24GeV^, the cross 
section is predicted to be aij^'N ^ ^2 + X)B{52 J/tpn^n^) = 4nb per nucleon. This 
value is in agreement with the preliminary result of E705, where the cross section is estimated 



to be 5.3 ± 1.9 ± l.Snb From Fig. 3, we can see that the experimental result of E705 
collaboration can not be explained within the color-singlet model, because the prediction of 
the color-singlet model is two orders of magnitude smaller than the color-octet model. This 
indicates that though the agreement with the E705 data can not be taken too seriously (due 
to large theoretical uncertainties, e.g., the choice of charm quark mass; the naive estimate for 
the color-octet matrix elements), an appreciable D-wave signal can only be explained by the 
color-octet mechanism. Experimentally, the strong signal of J /ipn^iT' at 3.836G'eV^ observed 
by i?705 is now questioned by other experiments. (E705 observed (77 ± 21) and (58 ± 21) 
background-subtracted events at ip'{2S) and 3.836 GeV respectively [14]; while E672-E706 
reported (224 ± 44) and (52 ± 30) background-subtracted events at ip'i^S) and 3.836 GeV 
respectively [20]). Nevertheless, if the £'705 result is confirmed (even with a smaller rate, 
say, by a factor of 3, for the signal at 3.836Gey, as might be implied by other experiments 
P0||), the color-octet production mechanism may provide a quite unique explanation for the 
D-wave charmonium production. 
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Moreover, the 1 D-wecve charmonium state ■0(3770) may also be observable at the 
fixed target experiments, with a rate smaller by a factor of | than the 2 state, provided 
the decay channel ■0(3770) — > DD (with about 100% decay branching ratio to the charmed 
meson pairs) can be detected. 

In conclusion, in this paper, we have calculated the D-wave heavy quarkonium production 
in fixed target experiments. We find that the color-octet mechanism plays an important role 
in the production. The color-octet contributions are two order of magnitude higher than 
the color-singlet contributions both in proton-nucleon and pion-nucleon collisions. Despite 
of some theoretical uncertainties, the prediction of the color-octet model is found to be in 
agreement with the preliminary result of E705 collaboration. This may provide another 
positive support to the color-octet production mechanism of heavy quarkonium. 
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APPENDIX: 



The short distance coefficients for the color-singlet processes gg ^Dj in Eq.( 

Figg ^ 'D^g) = '^""f^^^ — -{102M^'s'' + 302M'hH + 302M^%t^ 

+ 102M2°t^ - 286M^^s^ - 1732M^^s^t - 2844M^8s^t^ - 1732M^^st^ 

- 2mM^H* + 275MiS^ + 3840M^S^t + 10289M^S^t2 ^ i{)2mM^^sH^ 
+ 3840M^St^ + 275M^H^ - 227M^^s^ - 5004M^^s^t - 195Q9M^^sH'^ 

- 2953QM^^sH^ - 19569M^^sH'^ - 5004M^St^ - 227M^^t6 + AlOM^^s'^ 
+ 5137M^^sH + 23585M^^sh^ + mOSM^^sH^ + 47908M^^sh^ 

+ 23585M^2sV + 5137M^^st^ + 410M^V - 470M^°/ - 4220M^^sh 

- 1953AM^^sH^ - 47528M^°s^t^ - 63536M^°s^t^ - 47528M^°s^t^ 

- 1953AM^^sH'^ - 4220M^°st^ - 470M^°t^ + 245M^s^ + 2190M^s^t 
+ 10358M^sh^ + 28602M^sH^ + 47093M^sH'^ + 4:7093M^sH^ 

+ 28602M^sh^ + 10358M^sh'^ + 2190M^st^ + 2AbMh^ - 49MS^° 

- 580M*^s^t - 2822M^s¥ - 8984M^s¥ - 17653M^sV - 21968M^s¥ 

- l7Q^3M'^sH^ - 898AM^sH^ - 2822M^sh^ - 580MSt^ - 49M^t^° 
+ 67M^s^°t + 210M^s^t2 + 774M^s®t^ + 2006M^s^t^ + 3147M^s*'t^ 
+ 3147M^s^t^ + 2006MS^t^ + 774M^s3t^ + 210M^sH^ + 67MSt^° 
+ 25M^s^H^ + lOOM^sH^ + 220MhH'^ + 340MhH^ + 390M^sH^ 
+ 3A0M^sH'^ + 220M^sY + lOOM^sh^ + 25M^sY° + bs^^t^ 

+ 25sH^ + 60sh^ + 90sH^ + 90sH'^ + 60sh^ + 25sH^ + 5sh^^}, 

\yy 2y) 27 {M^ - sf{M^ -tf{s + tf^ 

- 3M^^/ - QAM^^sH - l3{}M^^sH^ - UM^^st^ - 3M^h'^ + 5M^S^ 
+ 180MiS^t + QllM^^sH^ + QllM^^sh^ + 180M^'^st^ - IIM^S^ 

- 28m^^sH - l3mM^^sH^ - 2208M^S^t^ - l3mM^^sh'^ - 28m^^st' 

- IIM^H^ + 20M^2^' + 269M^2^''^ ^ UlQM^'^sH'^ + 3927M^234^3 
+ 3927M^2^3^^ + 1716Mi2^¥ + 2691/12^^^ ^ 2{}M^H'' - 20M^°s^ 

- 1283M^°s^t2 - 3888Mi°s^t^ - 5450M^°sV - 3888M^°s^t^ 

- 1283M^^ sH^ - 144M^°st^ - 20Mi°t^ + lOM^s^ + IQM^sh + SM^^'t^ 
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+ 568M^sY + 16M^st^ + 10MH^2M^s^'^ + 20M^sH - UAM^'^sh 

- 156M^sh' - 1072M^sV - 2230M^s^t^ - 266AM^sH^ - 22S0M^sH^ 

- 1072MS^t^ - IbQM^sH^ + 20MSt^ - 2MH^^^ - 6M^s^°t 

+ SOM^s^t^ + A72M'^sh^ + 1172M^sV + 1570M^sh^ + 1570M^s^t^ 
+ 1172M^sH'^ + A72M^sH^ + bOM^sH^ - GM^st^^ - IGM^s^H' 

- leOM'sH^ - A96M'sH^ - S32M'sH^ - 960M'sH^ - 832M'sH'^ 

- 496M^s¥ - leOM'sH^ - IBM'sY'^ + 16s^H^ + SOsH"^ 
+ 192sH^ + 288sH^ + 288sH'^ + 192sH^80sH^ + 16sH^°}, 

+ 8M''H^ - 24M^^s^ - l28M^^sH - 206M^^sH' - 128M^^st^ 

- 2AM^H^ + 25M^S^ + SOOM^^sH + 826M^^sH'' + 826M^^sH^ 

+ SOOM^^st^ + 25M^H^ - UM^'Ks^ - SbGM^'^sh - IbbGM^'KsH'^ 

- 2A2AM^^sH^ - 1556M^^s2^^ _ 356M^^sf^ - UM^H!^' + im^'s^ 

+ 283Mi2si^ + lOM^^t^ - lOM^^s^ - l8QM^^sh - l2QlM^'^sH' 

- ?,U2M^'^sH^ - AmAM^'^sH'^ - ?,U2M^^sH^ - l2QlM^^sH^ 

- ISOM^^si^ - IQM^H^ + bM^s^ + 8QM^sH + 602M^sV 
+ 1943M^s^t3 + 3307M^s^t4 + 3307M^/t^ + l9A?,M^sH^ 

+ m2M^sh'^ + SOM^st^ + 5M^t^ - M^^^ - 20M^s9t - l98M^sh'' 

- 776M^sh^ - lb02M'^sH^ - 1812M^sH^ - lb02M^sH'^ - 776M^sH'^ 

- 198M^sH^ - 2{)M^st^ - MH'^^ + ?,M^s^H + A'^M^sH'^ + 22lM'^sH^ 
+ blAM'^sH'^ + 698MS^i^ + GQSM^s^i^ + blAM'^sH'^ + 221MS3i^ 

+ 3MSii° - SM^s^^i^ - SOM^s'^i^ - ISSM^s^i^ - 260M2s¥ 

- 300M=^/t^ - 260M2s^t^ - ISSM^s^t^ - SOM^s-V 

- hM'^sH^^ + Ss^^t^ + 25s9t4 + 60s®t^ + 90s¥ + 90s^t^ 
+ 60s^t^ + 25s¥ + 55^^^° + 283Mi2s^i}, 

where M = 2mQ and s, t, m are Mandelsterm invariants for the processes gg ^"Djg. 
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Figure Captions 



Fig.l. D-wavc charmonium production in proton- nucleon collisions for xp > 0. The dashed 
lines are the color-singlet contributions and the solid lines are color-octet contributions. For 
the solid lines, from up to down, they are for J = 3, 2, 1 D-wave sates respectively. For the 
dashed lines, from up to down, they are for J = 3, 1, 2 sates respectively. 

Fig. 2 D-wave bottomonium production in proton-nucleon collisions for xp > 0. The lines 
are the same as those in Fig.l. 

Fig. 3 D-wave charmonium production in pion-nucleon collisions for xp > 0. The lines are 
the same as those in Fig.l. 

Fig. 4 D-wave bottomonium production in pion-nucleon collisions for xp > 0. The lines are 
the same as those in Fig.l. 
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